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Abstract. We study the dependence of the intrinsic conductance of a nanocontact on its shape by using 
the recursion-transfer-matrix method. Hour-glass, torus, and spherical shapes are defined through analytic 
potentials, the latter two serving as rough models for ring-like and spherical molecules, respectively. The 
sensitivity of the conductance to geometric details is analyzed and discussed. Strong resonance efl'ects are 
found for a spherical contact weakly coupled to electron reservoirs. 

PACS. 73. 23. Ad Ballistic transport - 73.63.Rt Nanoscale contacts 



1 Introduction 



Nanocontacts are a research topic of intense current inter- 
est, motivated by the drive to reduce the size of integrated 
circuits Molecules, molecule-like species and nanopar- 
ticles have been considered as natural building blocks for 
interconnects. Es pec ially Benzene-dithiola te H, H , H, H, 0, 
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and 

SOjIMI have been extensively 



Ceo |23|,|M125| 
studied. 

Gaining detailed knowledge of bias-dependent quan- 
tum transport requires high-level ab initio calculations 
on the coupling of the molecular orbitals between the in- 
terconnect and the leads. However, an instructive view 
of ballistic properties of the contact can be achived via 
considering a general quantum mechanical transmission 
problem through an analytically defined constriction of 
various shapes. Here we calculated the conductance of a 
torus connected to jellium electrodes as a crude model 
for a Benzene-ring and the conductance of two concentric 
spheres which can be seen as a model for Cgo • In addition 
we investigated if two nanoconstrictions in close contact 
interact with each other. 

Among the variety of available methods to calculate 
quantum transport we have chosen the Recursion- Transfer- 
Matrix (RTM) method [H,!!!]. Compared to the Green's 
function methods which are widely used for calculating the 
conductance of nanocontacts the RTM method has the ad- 
vantage that it needs just a potential as input and uses 
no atomic orbitals. This feature made it possible to study 
simple potentials, only determined by geometry where pa- 
rameters can be changed easily. 



This paper is organized as follows: In Sect. [2] we briefly 
discuss the method. The results are discussed in Sect. [3] 
First we examine any interference effects when two nano- 
constrictions are in close contact. Then we investigate how 
the conductance of a torus depends on its geometry. Fi- 
nally we present results for the conductance of the concen- 
tric spheres for two cases: in the first case the constriction 
has cylindrical symmetry and in the second case the sym- 
metry is broken. Sect. S] concludes the paper. 



2 Method 

We calculated the conductance using the RTM method 
[HjllHl- In this method the Schrodinger equation is solved, 
where the potential is given on a grid and the energy E 
of the transmitting electron is chosen. As an ansatz for a 
constriction along the z-direction one expands the wave 
function into plane waves in the xy-plane 



E 



(t>tj{z) e 



iG\ 



(1) 



where G\ are reciprocal lattice vectors in the xy-plane 
and for the boundary condition one uses 



lentr 



(Z) 6''=^^+' 



(2) 

(3) 
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for extremal values of z, where (/)fj^^'^(z) denotes the wave 

function in the entrance and 4>^J^^^{z) denotes the wave 
function in the exit electrode. The components of the re- 
flexion matrix are denoted by r^. After the calculation 
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of the transmission matrix t 
Landauer-Biittiker formula 



— tij), one uses the 



(4) 



z=0: 



to calculate the conductance G, where G is measured in 
units of Go = 2e^ /h. The eigenchannels are the eigenval- 
ues of [111 . 

The convergence of the results was carefully tested 
with respect to the energy cut-off, the number of grid 
points and the distance between the grid points. We use 
atomic units throughout the paper. 



3 Results and Discussion 

3.1 Two constrictions in close contact 

To investigate if two constrictions in close contact interact 
with each other we calculated the conductance for the 
following geometry (see Fig. [T]): 

y = vb e{R + z) e{R-z) x 

X 0{x^ + {y+{R + W/2 + d/2)f - 

-(i? + W/2 - y/W^f) X 

X 0{x^ + {y-{R + W/2 + d/2)f - 



-{R + W/2^ ^R^ ~ z-^f) 



(5) 



We used Vo = 0.404, R^h imdW ^ 10. The distance 
between the constrictions was varied by varying d. We 
probed the contact at an energy E — 0.202. Our results 
show that there is no interaction between the two con- 
strictions. Additionally we made calculations using R = 2 
which are showing no interaction either. Due to the geo- 
metrical construction of the constriction it is not possible 
to calculate distances between the thinnest points of the 
constriction which are smaller than 2R + W . Anyway for 
smaller distances there may be chemical bonding if real 
atoms are considered. 



3.2 Torus connected to electrodes 

Another geometry which we studied is a torus. The ge- 
ometry we used is displayed in Fig. [2] where a cut along 
the yz-plane is shown. In addition the figure explains all 
the parameters. So c denotes the radius from the center 
of the ring to the center of the torus tube. The radius 
of the latter is called a. The torus is cut at z = ±D in 
a plane parallel to the xy-plane to open a pipe for the 
connection to the electrodes. This distance D also deter- 
mines the cross-section of the pipe since the pipe is just a 
prolongation of the hole which is created by the cut. The 
total length of this pipe is the sum of L and rp, where in 
addition rp is the inverse of the curvature of the opening 
from the pipe to the electrode. For these parameters we 
have chosen the following values: a = 4,c=10,£' = ll, 
L = 4, ro = 1, Fo = 0.404 and E = 0.202. 




Fig. 1. Geometry of the unit cell. The upper panel shows a cut 
through the contact perpendicular to the z-axis for z = 0. The 
lower panel shows a cut through the contact perpendicular to 
the a;-axis for a; = 0. At the white area the potential has the 
value V — Vo and the grey area denotes zero potential. The 
electrodes (in which V = 0) are located at z > R and z < —R. 




Fig. 2. Geometry of the unit cell at x = 0. At the white area 
the potential has the value V = Vo and the grey area denotes 
zero potential. The electrodes (in which V = 0) are located at 
z > D + L + ro and z < -~{D + L + ro). 
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Fig. 3. Conductance versus radius of the tube of the torus a. 
The other parameters have been chosen to be c = 10, D — 11, 
L = 4, ro = 1, Vo = 0.404 and E = 0.202. The insets show the 
geometry of the contact perpendicular to the a;-axis for a = 2 
(left inset) and a = 10 (right inset). 



Fig. 4. Conductance versus D, where _D is a measure not 
only where the torus is cut but also for the size of the pipes 
connecting the torus to the electrodes. The insets show the 
geometry of the contact perpendicular to the a:;-axis for D — 6 
(left inset) and D = 13 (right inset). The other parameters 
have been chosen to be a = 4, c = 10, L — 4, ro — 1, Vo ~ 0.404 
and E = 0.202. 



To investigate the dependence of the conductance on 
these geometrical parameters, calculations have been done, 
where we varied one parameter and left the other param- 
eters constant: Figure [3] displays conductance against the 
radius a of the tube of the torus. The larger a the larger is 
the conductance. Increasing a increases also the diameter 
of the pipes, which connect the torus to the electrodes, 
since D is kept constant. 

One can see conductance steps ata = 3,a = 5.5 and 
a = 7.5. This can be understood as follows: In the calcula- 
tion an electronic energy of 0.202 is used. This corresponds 
to a wavelength A of 9.9. When a is smaller than about 
2.7 then the tube and the pipes are too narrow to sup- 
port half a wavelength and thus there is no current going 
through. For 2.7 < a < 5.1, half a wavelength fits through 
the constriction and thus we can see the first conductance 
step. If 5.1 < a < 7.4, the system can support one wave- 
length and the conductance jumps up to G = 5 Gq. The 
next conductance step takes place at a = 7.4 because from 
a = 7.4 on until a = 9.9 the constriction is able to support 
1.5 wavelengths. If a is larger than 9.9 two wavelengths are 
fitting inside the system and so one would expect the next 
conductance step at a = 10. 

The dependence of the conductance on D is displayed 
in Fig. m As expected the conductance decreases with a 
decreasing cross-section of the pipe (increasing D). 

Fig. [5] displays the energy-dependence of the conduc- 
tance for the considered geometry. For energies larger than 
about 0.1 the system conducts. The increase of the con- 
ductance with increasing energy goes stepwise. The second 
conductance step takes place around E = 0.3. At the con- 
ductance plateaus the conductance shows oscillations. At 
the first conductance plateau the conductance oscillates 
between 1 and 2 Go, on the second plateau the conduc- 
tance takes values between 4 and 5 Gq. 
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Fig. 5. Conductance versus energy E for a geometry where 
long pipes connecting the torus to the electrodes. The inset 
shows the geometry of the contact perpendicular to the x-axis. 
For this calculation the following parameters have been chosen: 
a = 4, c = 10, D = 11, L = 4, ro = 1 and Vo = 0.404. 



To understand the conductance steps one can argue in 
a similar way than when varying a (Fig. [3]). In this case a 
is fixed to 4, which means that the wavelength has to be 
shorter than 16 (8) in order to make half a (one) wave- 
length fit inside the tube. The connection to the electrodes 
has an extension of 7.75 in x-direction and an extension 
of 16.61 in y-direction. Converting these lengths in en- 
ergy yields that the system should start to be conducting 
around E = 0.1 and the next conductance channel should 
open around E = 0.3. 

The oscillations are an effect of the long pipes between 
the electrodes and the torus. Refs. [sH] and [S^ calculate 
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Fig. 6. Conductance versus energy E for a geometry where 
the torus is almost directly connected to the electrodes. The 
inset shows the geometry of the contact perpendicular to the 
X-axis. For this calculation the following parameters have been 
chosen: a = 4, c = 9, D = 10, L = 0, ro = 0.5 and Vo = 0.404. 



narrow pipes and find similar oscillations the longer the 
pipes. 

Figure [6] shows conductance versus energy for another 
set of parameters. Here the parameters are a — A, c — 9, 
D = 10, L = 0, ro = 0.5 and Vq = 0.404. This means 
that now the total length of each pipe is one order of 
magnitude smaller than in Fig. [5l The differences in the 
other parameters are small, which means that the main 
difference is due to the length of the pipes which connect 
the torus to the electrodes. 

We observe that the conductance steps take place at 
the same energies than in Fig. [51 This is due to the same 
size of the tubes. In both cases a = 4 and the exten- 
sion of the pipe connecting the torus to the electrodes in 
x-direction is in both cases 7.75. Only the extension of 
these pipes in y-direction is here broader (17.32 instead 
of 16.61). But this time there are much less oscillations at 
the plateau with a much smaller amplitude. This supports 
the theory that the oscillations may be due to the length 
of the straight pipes. 

Though our calculations are only a coarse model, we 
can nevertheless compare them with experiment. There 
the conductance of a benzene-l,4,-dithiolate molecule be- 
tween gold electrodes has been measured [1^1 and a step- 
like structure of G(V) has been found. Since G{E) can 
be seen as an approximation to G{V), our results agree 
qualitatively with experiment. 



3.3 Sphere connected to electrodes 

The geometry used for calculating the conductance of a 
sphere is displayed in Fig. [T] Due to the symmetry cylin- 
drical coordinates have been used where r = \J + y^- 
At the white area the potential has the value V = Vq 
and the grey area denotes zero potential. Two concentric 
spheres with radii Ri (inner sphere) and Ra (outer sphere) 





Fig. 7. The geometry used to calculate the conductance 
through a sphere connected to electrodes. Due to the symme- 
try cylindrical coordinates have been used. At the white area 
the potential has the value V = Vo and the grey area denotes 
zero potential. The electrodes (in which V = 0) are located at 
z > D + L + ro and z < -{D + L + ro). 



have been used in order to be able to change the thickness 
of the spherical layer where the electrons are allowed to 
flow. At a distance D from the center of the spheres in ±z- 
direction the sphere is cut parallel to the xy-plane. This 
cut determines the radius of the pipes which connects the 
sphere to the upper and lower electrode via 



(6) 



The total length of each pipe is L + r^ where the inverse of 
To is denoting the curvature of the opening from the pipe 
to the electrode. 

The first panel of Fig.[S]shows the conductance plotted 
against energy for a sphere which is weakly coupled to the 
electrodes. In this case "weak coupling" means that there 
is only a small connection between the sphere and the 
electrodes. The parameters used in this plot are Ri = 1.7 
and Ra = 7.7 for the radii of the concentric spheres, D ~ 
7.6, L — 0.0 and ro = 0.5. Thus there is only a small hole 
of radius ro = 1.2 connecting the sphere to the electrodes. 

For conduction this small hole must be able to contain 
at least half a wavelength. Using 



E 



2^ 
A2 



(7) 



for conversion between wavelength and energy, the system 
should start to conduct at 0.86. This energy is out of the 
range displayed in Fig. [51 But what is seen in Fig. [51 are 
resonances at special energies. These resonances display 
the eigenstates of the sphere. 
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Fig. 8. Top panel: Conductance versus energy for a sphere 
which is weakly coupled to the electrodes. The three lower 
panels show the first three eigenchannels. A cut of the potential 
along the yz-plane is shown in the inset. The parameters used 
are Ri = 1.7 and Ra = 7.7, D = 7.6, L = 0.0 and ro = 0.5. 



The three lower panels of Fig. [8] show the first three 
eigenchannels of the eigenchannel analysis for the geome- 
try used in Fig. [5] The first channel is non-degenerate but 
the second and third eigenchannels are degenerate due to 
the cylindrical symmetry of the constriction. 

A different geometry has been used to calculate the 
conductance versus energy curve shown in the first panel 
of Fig.m Here Ri = 3.7, Ra = 9.7, D = 8.6, L ^ 0.0 and 
ro = 1.0. Now td = 4.5 and using the same argumenta- 
tion than above one expects that the system conducts for 
energies larger than 0.06. Fig. [9] confirms our prediction: 
around E = 0.06 the conductance jumps from Gq to 
1 Gq. With increasing energy the conductance increases 
but sharp peaks are superpositioned on this increasing 
conductance. These sharp peaks again display the eigen- 
states of the sphere. An eigenchannel analysis has been 
made and the first three eigenchannels are displayed in 
the lower panels of Fig. [51 The first eigenchannel opens 
at 0.06 and stays open for higher energies showing some 
oscillations. The second and the third one are degenerate 
and are showing some resonances before they open around 
an energy of 0.2. 




E / Hartree 

Fig. 9. Top panel: Conductance versus energy for a sphere 
which is strongly coupled to the electrodes. The three lower 
panels show the first three eigenchannels. A cut of the po- 
tential along the yz-plane is shown in the inset. Geometrical 
parameters used are Ri = 3.7, Ra = 9.7, D = 8.6, L = 0.0 and 
ro = 1.0. 



The geometries used for obtaining the results in Fig. 
[8] and Fig. [9] have almost no pipe, the spheres touch more 
or less directly the electrodes. What changes if the con- 
nection to the electrodes is made by long pipes? Fig. [TO] 
shows results obtained for a geometry where all param- 
eters are the same than in Fig. [9] except L. This time 
L = 5.0 is used instead of L = 0.0. Now the picture looks 
much more complex. G{E) (first panel) shows resonances 
and antiresonances superposing a "background" conduc- 
tance. The complexity may arise due to the much more 
complex geometry since every change of geometry (from 
electrode to pipe, from pipe to sphere and so on) gives 
rise to reflections. The three lower panels of Fig. [TO] show 
the first three eigenchannels of which two are degenerate. 
This time the second and third eigenchannel open at an 
energy of about 0.24. The opening of these two channels 
can be seen in the first panel of Fig. [TO] when the average 
conductance rises above one conductance quantum. 

One can make a qualitative comparison to experiments 
where the conductance of a Cgo-molecule has been mea- 
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Fig. 10. Top panel: Conductance versus energy for a sphere 
which is strongly coupled through long pipes to the electrodes. 
The three lower panels show the first three eigenchannels. A 
cut of the potential along the yz-plane is shown in the inset. 
Geometrical parameters used are Ri = 3.7, Ra = 9.7, D = 8.6, 
L = 5.0 and ro = 1.0. 

sured [30'|. In these measurements a peaked structure of 
the conductance has been found as well. 



3.4 Sphere connected asymmetrically to the electrodes 

Until now cylindrical symmetry has been assumed for the 
constriction. In addition all constrictions had a mirror 
plane in the xy-plane. How does it affect the conductance 
if these symmetries are broken? 

To answer these questions calculations have been done 
where one pipe was shifted away from the z-axis, where yp 
denotes the distance from the 2;-axis and Xp — without 
loss of generality. Fig. [11] shows a cut of the geometry 
along the yz-plane. 

Fig. [12] shows conductance versus yp. The other pa- 
rameters were chosen to be i?a — 10, Ri = 5, P — 11, 
L = 2, Z? = 8, rp = 3, ro = 2, Vb = 0.404, E = 0.202. 
If the upper pipe is situated at the z-axis (yp = 0), the 
conductance is about 0.82 Gq. If the upper pipe is shifted 
away from the z-axis, the conductance increases until it 
reaches a maximum at yp ~ 2. There the conductance is 



through analytic constrictions 
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y 



Fig. 11. The geometry used to calculate the conductance 
through a sphere connected asymmetrically to electrodes. At 
the white area the potential has the value V = Vq and the grey 
area denotes zero potential. The electrodes (in which V = 0) 
are located at 2: > P -I- ro and z < —{D + L + ro). 




Fig. 12. Conductance versus yp. The other parameters take 
the following values: Ra = 10, Ri = 5, P = 11, L = 2, D = 8, 
rp = 3, ro = 2, 1/0 = 0.404, E = 0.202. 

0.84 Gq. ?/p = 2 is special in that way, that at this ge- 
ometry the outer edge of the pipe has the same value of 
the y-coordinate than the surface of the inner sphere. If 
the distance of the upper pipe to the z-axis is increased 
further, the conductance decreases. 



4 Conclusions 

We calculated the conductance of nanocontacts with dif- 
ferent geometries using the recursion-transfer-matrix method. 
The calculated geometries are two constrictions in close 
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contact, a torus and two concentric spheres. Our results 
show that geometry matters. It is not only the diameter of 
the smallest part of the constriction, which determines the 
conductance. Furthermore small changes in the geometry 
can have a large impact on conductance. 

For the case of two constrictions in close contact we 
found, that they do not interact with each other. 

In the case of the torus we found a steplike behaviour of 
the conductance. Each time when there fits half a wave- 
length more inside the tube by either increasing the ra- 
dius of the tube or the energy of the electrons there is a 
conductance-step. 

For the sphere a discrete set of energies - its eigenstates 
- is obtained if the coupling to the electrodes is weak. In 
the case of strong coupling the resonances remain but the 
overall behaviour resembles more a step structure. 

Moreover the way how the torus or the tube is con- 
nected to the electrodes is important for the conductance. 
Connecting the torus or the sphere with pipes to the elec- 
trodes causes more fluctuations in the conductance since 
the electrons are suffering more reflections. The total con- 
ductance is decreased compared to the case with pipes of 
negligible length. 

We would like to thank Matti Manninen for fruitful discussions 
and the Academy of Finland for financial support. 
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